JV/zA/sf ' 

/p - / 8> 


75 <5 / 


f y *fl 

A Robust Momentum Management 

AND A TTITUDE CONTROL SYSTEM FOR 

the Space Station 


J.L. Speyer 
Ihnseok Rhee 


University of California , Los Angeles 


imiSSe J for 

THE SPACE STATION (Houston Omv.) ^ 


N9I-20202 


Unci ds 
0007301 


Cooperative Agreement NCC 9-16 
Res earch Activ ity No. MS.4 


NASA Johnson Space Center 
" Engineering Directorate 

- Navigation, Control, & Aeronautics Division 



Research Institute for Computing and Information Systems 
University of Houston - Clear Lake 



TECHNICAL REPORT 



The 

RICIS 

Concept 


The University of Houston-Clear Lake established the Research Institute for 
Computing and Information systems in 1986 to encourage NASA Johnson Space 
Center and local industry to actively support research in the computing and 
information sciences. As part of this endeavor, UH-Clear Lake proposed a 
partnership with JSC to jointly define and manage an integrated program of research 
in advanced data processing technology needed for JSC’s main missions, including 
administrative, engineering and science responsibilities. JSC agreed and entered into 
a three-year cooperative agreement with UH-Clear Lake beginning in May, 1986, to 
jointly plan and execute such research through RICIS. Additionally, under 
Cooperative Agreement NCC 9- ^computing and educational facilities are shared 
by the two institutions to conduct the research. — — 

The mission of RICIS is to conduct, coordinate and disseminate research on 
computing and information systems among researchers, sponsors and users from 
UH-Clear Lake, NASA/JSC, and other research organizations. Within UH-Clear 
Lake, the mission is being implemented through interdisciplinary involvement of 
faculty and students from each of the four schools: Business, Education, Human 
Sciences and Humanities, and Natural and Applied Sciences. 

Other research organizations are involved via the “gateway” concept. UH-Clear 
Lake establishes relationships with other universities and research organizations, 
having common research interests, to provide additional sources of expertise to 
conduct needed research. 

A major role of RICIS is to find the bes t m atch of sponsors, researchers and 
research objectives to advance knowledge in the computing and information 
sciences. Working jointly with NASA/JSC, RICIS advises on research needs, 
recommends principals for conducting the research, provides technical and 
administrative support to coordinate the research, and integrates technical results 
into the cooperative goals of UH-Clear Lake and NAS A/ JSC. 


A Robust Momentum Management 

AND A TTITUDE CONTROL SYSTEM FOR 

the Space Station 







Preface 


This research was conducted under the auspices of the Research Institute for 
Computing and Information Systems by J.L. Speyer and Ihnseok Rhee of the 
University of California, Los Angeles. Dr. A. Glen Houston, Director of RICIS, 
served as RICIS research representative. 

Funding has been provided by Navigation Control & Aeronautics Division, 
Engineering Directorate, NASA/JSC through Cooperative Agreement NCC 9-16 
between NASA Johnson Space Center and the University of Houston-Clear Lake. 
The NASA technical monitor for this activity was David Geller, of the Navigation 
Section, Navigation and Gridance Systems Branch, Navigation Control & 
Aeronautics Division, Engineering Directorate, NASA/JSC. 

Hie views and conclusions contained in this report are those of the author and 
should not be interpreted as representative of the official policies, either express or 
implied, of NASA or the United States Government. 





A ROBUST MOMENTUM MANAGEMENT AND 
ATTITUDE CONTROL SYSTEM FOR THE SPACE 

STATION 


rhnseok Rhee* and Jason L. Speyer^ 

The University of Texas, Austin, Texas 78712 


Abstract 

A game theoretic controller is synthesized for momentum management and atti- 
tude control of the space station in the presence of uncertainties in the moments of in- 
ertia. Full state information is assumed since attitude and attitude rates are assumed 
to be very accurately measured. By an input-ouput decomposition of the uncertainty 
in the system matrices, the parameter uncertainties in the dynamic system are repre- 
sented as an unknown gain associated with an internal feedback loop(IFL). The input 
and output matrices associated with the IFL form directions through which the un- 
certain parameters affect system response. If the quadratic form of the IFL output 
augments the cost criterion, then enhanced parameter robustness is anticipated. By 
considering the input and the input disturbance from the IFL as two noncoopera- 
tive players, a linear-quadratic differential game is constructed. The solution in the 
form of a linear controller is used for synthesis. Inclusion of the external disturbance 
torques results in a dynamic feedback controller which consists of conventional PID 
control and cyclic disturbance rejection filters. It is shown that the game theoretic 
design allows large variations in the inertias in directions of importance. 
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I. Introduction 


A game theoretic controller developed in Ref. 1 is applied to the attitude/momentum 
control for the space station which uses control moment gyros(CMGs) as the primary 
actuating devices and gravity gradient torque to manage momentum stored in CMGs. 
The moments of inertia of the space station are assumed constant but uncertain. In 
Ref. 2,3 the linear quadratic regulator(LQR) design procedure has been used to con- 
trol the attitude/ momentum of the space station. Full state information is assumed 
since the attitude and attitude rate are assumed to be very accurately measured. In 
Ref. 2 disturbance rejection filters are augmented to the system to handle the exter- 
nal cyclic disturbance torque, and the LQR design and pole assignment procedures 
for pitch control and roll-yaw control, respectively, are applied to the augmented sys- 
tem. In this paper the system equation is differentiated until the external disturbance 
torque term disappears in the resulting equation in order to apply the design proce- 
dure developed in section II. The resulting controller consists of conventional PID 
control and the cyclic disturbance rejection filter as in Ref. 2. 

The application of the game theoretic approach combined with the internal feed- 
back loop decomposition for describing parameter uncertainty allows very large vari- 
ation in the inertia of the space station with little deterioration in performance. In 
Ref. 4, a differential game approach to developing synthesis techniques was taken 
where the parameter uncertainty was not decomposed and only the uncertainty in 
the system matrix is considered. In Ref. 5-7, Lyapunov stability theory has been 
used to design a control law for a system with uncertainty. This approach is similar 
to that used here in a particular algebraic Riccati equation(RDE). In Ref. 8, by ad- 
dopting an input-output decomposition of the parameter uncertainty, the uncertain 
system is represented as an internal feedback loop(IFL) in which the parameter un- 
certainty is embedded in the system as a fictitious disturbance. Tahk and Speyer 8 
developed the parameter robust linear-quadratic Gaussian(PELQG) synthesis pro- 
cedure which is an LQG design based on an extension of loop transfer recovery for 
the IFL description. In Ref. 7,8 the system is augmented to accommodate the input 
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matrix uncertainty. Theoretically, this approach is limited in that input and output 
matrices associated with the IFL are to have the dimension of the orignal input and 
outputs, respectively. In the game theoretic approach, this restriction is not required. 
By considering the input and fictitious input in the IFL description as two nonco- 
operative players, a finite-time linear differential game problem is constructed. By 
taking the quadratic norm of the fictitious output, the cost criterion is augmented by 
a term which emphasized robust performance. By taking the limit to an infinite- time, 
time- invariant linear system, a time- invariant control law is obtained. It is shown that 
the resulting time- invariant controller stabilizes the uncertain system for a prescribed 
parameter uncertainty bound. The development of the game theoretic controller is 
presented in section II. The approach taken in Ref. 1 generalizes the results here to 
the partial information problem where only some noisy measurements of the states 
are available. 

One motivation for this paper is to demonstrate on a meaningful problem the 
design process using the game theoretic controller augmented with the IFL decompo- 
sition. Although the linear-quadratic regulator has guaranteed gain an phase margin, 
many systems remain sensitive to parameter variations. This control problem is par- 
ticularly interesting in that the variation in the moments of inertia are bounded by 
physical constraint. The IFL decomposition allows selective changes in the moments 
of inertia to be included in the design process. In the pitch channel, there are two 
independent parameter uncertainties, one associated with the system matrix and the 
other associated with the input matrix. The quadratic norm of the fictitious output 
from the IFL decomposition of the system matrix augments the quadratic cost crite- 
rion and this augmented term represents a measure of system robustness. The effect 
of the fictitious inputs by the decomposition of the input matrix is to increase the gain. 
However, in the roll-yaw axis where there are three independent parameters, stability 
robustness in directions associated with inertia variations that -an be made large 
before reaching physical constraint is achieved without increased bandwidth. The es- 
sential design task is choosing the weighting for combining the parameter uncertainty 
directions which improve stability robustness subject to the physical constraints on 

3 





the inertias. 

This work is based upon reports Ref. 9,10 where the internal feedback loop 
concept 8 was applied to the problem of improving robustness of momentum manage- 
ment and attitude control in the roll-yaw axes using LQR theory. The game theoretic 
controller 1 first suggested in Ref. 11 was first applied to this problem in Ref. 12 for 
the pitch axis. At that time the authors became aware of the work in Ref. 13 for the 
roll-yaw axes using similar techniques, and latter for all axes in Ref. 14,15. 

II. Game-theoretic controller 

A controller for a linear time-invariant system with parameter uncertainties in 
the system and input matrices is derived via the differential game frame work. A 
game theoretic approach is taken because it is shown under certain conditions that 
there exist a nonnegative definite solution to an ARE, then the disturbance atten- 
uation function is bounded 1 . This is equivalent to imposing an norm bound 
on the transfer function between the disturbance input and the desired output. 1 In 
this section the disturbance inputs associated with system parameter uncertainty are 
constructed by the internal feedback loop decomposition of Ref. 7,8. 

Consider a time-invariant linear system with uncertainties in the system and input 
matrices described by 

* = (Ao + AA)x + (B 0 + AB)u ( 1 ) 

where x € R n , u € RT, A 0 G #**", and B 0 6 S" Xm denote the state, the input, 
the nominal system matrix, and the nominal input matrix, respectively, and AA and 
AS are perturbations of the system matrix and the input matrix, respectively, due to 
parameter variations. It is assumed that all states are directly measured and (Aq, So) 
is a stabilizable pair. 

By adopting the input-output decomposition modeling 3 of the perturbations, A A 
and AS, are represented as 

A A = DL a (e)E, A B = FL b (e)G (2) 
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where e denotes the parameter variation vector which is constant but unknown, and 
D, E, F, and G are known constant matrices. It is noted that the elements of £ 
need not be independent of e; ch other. With this modeling of A,4 and AS, the 
uncertain dynamic system (1) can be represented as an internal-feedback-loop(IFL) 
description 8 in which the system is assumed forced by fictitious disturbances caused 
by the parameter uncertainty: 


x = ,4 0 x + B 0 u + T w 


yi = 


x + 


0 

G 


where w = 


E 

0 

wj _ L a (e) 0 

w 2 0 L b {e) 

[ w i" is the fictitious disturbance and r = [D /*]. 


u 


yi 


( 3 ) 

( 4 ) 

( 5 ) 


In the above IFL description the fictitious disturbance, w, is a feedback signal of 

r L a (e) 0 
0 L b (e) 


yi amplified by the unknown gain 


. Hence, one way to reduce the 


effect of parameter uncertainty is to assume that w is an independent Gaussian white 
noise and design a controller minimizing the cost 


.I'Sej; l’ il>yTy+y ' y '' >dt 

subject to the system equation (3) where y is a performance measure defined as 



c 


0 

y = 


x + 



0 J 


. Cl . 


and p is a positive constant which represents the trade-off between the performance 
described by y T y and the robustness with respect to parameter uncertainty described 
by yfyi- Let 


Q = pC T C + E t E 
R = pCf Ci + G t G. 

By assuming ( Q , ,4o) is detectable, R is positive definte, and G = 0, this cost criterion 
leads to the PRLQG design procedure 8 as p — ► 0. 
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An alternate approach to robust synthesis is to design a controller to make the 



disturbance attenuation function due to the fictitious disturbance bounded, 

sup 

u/e^[o.^] Jq w 1 w at 

i.e., 






where 7 is a positive constant, and tj is a fixed final time. This problem can be solved 
by solving a differential game 16 to find u that minimizes and w that maximizes the 



cost criterion 


— 


(y r y + yfyi - 7 _ 2 w r w) 

(7) 

m 


subject to (3). It is well-known 4,17 that if there exists a real symmetric solution FI(£) 

W 


over the interval t G [0,*/] to the Riccati differential equation(RDE) 




-n = + nAo - n (b 0 r- 1 b* - 7 2 rr r )n + q 




with the final condition II(£/) = 0 , then the strategies for u and w described 

as 



U*= -R~ 1 B T l 1(0 X 


— 


w* = 7 _ 2 r r n(t)x 




yield the saddle point, i.e., 


III 1 


J(u m ,w,tf) < ) < J(u,w m ,t f ) 

( 8 ) 

— 


Vu,w 6 L 2 [0,t f ]. 


w 


For the case where tj — ► oo, II(i) converges to a constant matrix if there exists a 

— 


nonnegative definite solution to the algebraic Riccati equation(ARE) 




o = Alti + fUo - fUBoR-'B? - 7 2 rr T )fi + q. 

(9) 

Ill 


Note that in general there may be many nonnegative definite solutions to the ARE 
(9). The minimal nonnegative definite solution 1,6 to the ARE (9), denoted as IT, is 

Ip R 


defined as a nonnegative definite solution to the ARE (9) such that fl < ft where fl 
is any nonnegative definite solution to the ARE (9). Then Il(<) — ► ft as tj -+ co . 1,17 



Hence, u* and w* become time- invariant strategies described by 


III II 


u= -R~ 1 B T Ux 

( 10 a) 



w = 7 _ 2 r T n x 

( 10 b) 
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The resulting time-invariant strategies (10), however, may not satisfy the right hand 
inequalities in (S) as tj — * oo. 1. However, only the left hand inequality is of concern 
in the development of this class of robust controllers. 

In the worst case design, since the fictitous disturbance w is not an intelligent 
player, only the control strategy for the control u given by (10a) can be implemented. 
The following proposition provides a robustness property for the control law (10a). 

Proposition 1 Assume that R is a positive definite matrix and (Q,A 0 ) is a de- 
tectable pair. Suppose that there exists a nonnegative definite solution, fi, to the 
ARE (9). Then, the control law given as 

u = - R~ l B T fix (11) 

stabilizes the uncertain dynamic system (1) for all e such that ||Z, a (e)|| < 7, and 
11^)11 < 7- 

Claim 1 Suppose that Dj D x + G T G > 0. Then, 

DfU l D 1 +G T U 2 G > 0 VU U U 2 > 0. 

Proof) It is sufficient to prove DjU^Di + G T U 2 G is nonsingular. Suppose that there 
exists a nonzero vector, z, such that 

z T (DjU l D l + G T U 2 G) z = 0. 

Then, D\Z = 0 and Gz = 0 since U\ and U 2 are positive definite, hence (Dj D\ + 
G T G)z = 0 which contradicts the assumption. □ 

Proof of Proposition 1) By using the control law (11), the closed loop system is 
discribed as 

x = A c x (12) 

where 


A c = A 0 + DL a (e)E - {B 0 + FL b (£)G)R- l B^tl. 
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lAOS’-S'fl' 



The ARE (9) can be rewritten as following the Lyapunov equation: 

AlYl + UA C = - Q x (13) 

where 


Q i = nB Q R- x k h R- x Bin + £ r A a £ + pC T C 
+ 7 2 (nz? - r 2 E T Ll)(tiD - 1 ~' 1 e t lI) t 
+ 7 _2 n( 7 2 F + B 0 R-'G T Ll)(~f 2 F + B 0 R- 1 G T Ll) T n 
A a = I-r 2 L a (e) T L a (e) 

A b = P ClC 1 + G t (I - 7 ~ 2 L b (e) T L b (e))G. 

||L 0 (e)|| < 7 implies that A a > 0, and ||X 6 (e)|| < 7 and claim 1 yield A 6 > 0. Hence, 
Qi is nonnegative definite. Now it will be shown that (Qi, A e ) is a detectable pair by 
contradiction. Suppose (Qj, A c ) is not detectable. Then, there exists a nonzero vector 
z for some s in the closed right half plane such that (si — A c )z = 0 and Q x z = 0. 
Since each term in Q x is nonnegative definite, z T Q x z = 0 leads to 

z T (tlB Q R- 1 A b R- 1 Bln + E T A a E + pC T C) z = 0 


which implies that B% tlz = 0, Ez = 0, and Cz = 0, hence 


(si — A c )z = (si — Ao)z. 


Therefore, 


si — Ao 
pC T C + E T E 


z = 0 


which contradicts to the assumption (Q,Ao) detectable. Applying the lemma 4.2 18 
to the Lyapunov equation (13) completes the proof. □ 

Note that proposition 1 holds for any nonnegative solution to the ARE (9). However, 
the minimal nonnegative solution, IT, produces the smaller gain for the control law. 
In order to design the controller (11), the design parameters p and 7 should be 
chosen for the ARE (9) to have a nonnegative definite solution. In particular, as 
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the value of p increases, system performance improves but the stability robustness 
with respect to the parameter variation becomes poor. As the value of 7 increases, 
stability robustness with respect to parameter variation improves. 


■'w 


III. Space Station Control 


The game theoretic controller developed in section II is applied to the attitude/momentum 
control for the space station. 

A. Space Station Dynamics 

The space station is expected to maintain a local vertical/local horizontal(LVLH) 
orientation during normal operation. Suppose that the space station control(body) 
axes are aligned with the principal axes. (For the phase-I configulation of space station, 
this is a good assumption. 2 ) For the small deviation from LVLH frame, the linearized 
space station dynamics are described as 2,3 


^ + 4 w 2 A: r <£ - iu 0 (l - k x )ip = - y{T x - w x ) 
6 - Zu]k y 0 = — y{T y — Wy) 
V* - + uj 0 ( 1 + k z )<j> = — j ~{T Z — w z ) 

* Z 


(14a) 

( 14 b) 

( 14 c) 


where the body fixed axes (1, j/,2) denote the roll, pitch, and yaw control axes with 
the roll axis in flight direction, the pitch axis normal to the orbit plane, and the yaw 
axis toward the Earth; <j>, 0, and tp denote the roll, pitch, and yaw Euler angles with 
respect to the LVLH frame; (T x ,T y ,T z ) is the control torque vector produced by the 
CMG with respect to the control axes; (w x , w y ,w z ) is an external disturbance torque 
vector with respect to the control axes; uj 0 is the orbital rate of 0.0011 rad/sec; and 
k x , k y , and k z are the parameters defined from the moments of inertia, I x , I y , and I z 



Terms involving w 2 in ( 14 ) represent the combined gravity gradient and gyroscopic 
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torque in each axis. The CMG momentum dynamics are' 


(16a) 

(16b) 

(16c) 


2,3 

h x h. - — Tj- 

hy — Ty 

h z + u 0 h x — T, 

where (h x , h y , h x ) are the CMG momentum vector with respect to the control axes. 
It is assumed that the Euler angle, the Euler angle rate, and the CMG momentum 
are perfectly measured. The roll and yaw dynamics are coupled while the pitch axis 
is uncoupled. 

The physical constraints for the parameters k x , k y , and k z due to the triangular 
inequality of moment of inertia 19 are 

i = z,y,z. (17) 

The moments of inertia, 7 X , I v , and I z are assumed constant but uncertain and de- 
scribed by 

Ii=kn + A/,-, i = x,y,z (18) 

where the subscript ‘n’ and A Ii denote the nominal value and the variation of each 
moment of inertia, respectively. Then, the parameters, k x , Ar y , and k z can be repre- 
sented as 


ki = k in + A k,, i = x, y, z (19) 

where denotes the value of ki with nominal values of the moments of inertia and 
A A:,- denotes the variation due to the variation of the moments of inertia. The nominal 
values of the moments of inertia for the Phase- 1 configuration are 

I xn = 50.2SE6, I yn = 10.S0E6, I 2n = 58-57E6 

in unit of slug-ft 2 . In order to check a stability margin for inertia variation, ten types 
of variations of the moments of inertia listed in Table 1 are considered. Each variation 
is limited by the physical constraint (17). Table 1 also shows the physical limit of 
each type of variation. 
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The external disturbances (w z ,w y ,w z ) are modeled as 2,3 

w i = A d u sin (u 0 t + ip u ) + /if, sin(2u> 0 t + <^ 2> ) + Bf, (20) 

i = x,y,z where /if,, /if,, and Bf are assumed constant but unknown. The cyclic 
aerodynamic disturbance at orbital rate and twice the orbital rate are due to the 
diurnal bulge and the rotating solar panel, respectively. 

B. Pitch Control 

Before developing the controller for the pitch-axis, the open loop characteristics 
of the pitch channel are investigated. For the external disturbance free case, suppose 
the constant feedback control described as 

9 

9 

h y 

stabilizes the nominal system of (14b) and (16b). By using this control the closed 
loop characteristic equation, A(s), becomes 

A(s) = s*-(k 3 - jz'j s 2 + (^- -3u 2 0 k y y+3u 2 0 K 3 k y . 

By the Routh-Hurwitz criterion 21 , the zeroth-order term in right hand side should 
be positive for k y = k yn . Hence, the given feedback control law can not stabilize the 
system when the sign of k y is different from that of k yn . In other words, any constant 
feedback control law designed for the nominal system can not stabilize the system 
with a k y whose nominal value has a different sign. 

The appearance of a cyclic disturbance, described in (20), prevents the direct 
application of the game theoretic controller to the pitch channel. However, this 
can be avoided by differentiating (14b) and (16b) until the cyclic disturbance term 
disappears in the resulting equation. Differentiating (14b) and (16b) five times yields 



9 - (3 k y - 5 )u#< v > + (15 k y - 4 K0< 7// > 

(21) 


+ 12 kyiJ^O — fyUy 


9 

h { y VI) = - bu]h[ IV) - 4 

(22) 
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where the parenthetical superscripts represent the order of the time derivative, f y — 
If-, and Uy is a new control variable defined as 

; y 

“v = T- {T< V) + 5 + 4u ,tt y ) (23) 

i yn 7 

Note that the parameter f v has uncertainty due to the uncertainty in I y and is rep- 
resented as 

fy = 1 + A/y (24) 

Equations (21) and (22), however, are not yet an adequate representation of the 
system equation for the design procedure developed in section II since they contain 
uncontrollable modes on the imaginary axis when u y is used as a control. It can be 
verified that the uncontrollable modes are at s = 0, s = ±ju> 0 and s = ±2 ju 0 which 
arise from differentiation of the cyclic disturbance. The uncontrollability problem can 
be avoided by changing the regulated variables. The uncontrollable mode at s = 0 
can be removed by regulating 9 instead of 9. If h y is regulated instead of h y , h y 
becomes unbounded as time increases. It is clear from original pitch dynamics and 
CMG momentum equations that the pitch attitude 9 can not be regulated since T y in 
(14b) requires a biased control to regulate 9 in steady state under the disturbance^ujy. 
Hence, h y becomes unbounded. Note that regulating h y instead of h y still produces 
the uncontrollable mode at s = 0. Since the uncontrollable oscillating modes at 
s = ±ju 0 and s = ±2ju 0 arise in ( 21 ) and ( 22 ), all the modes of the 9 and h y 
channels can not be regulated. However, as will be shown, the oscillating modes in 
either the 9 or h y channel can be regulated. To regulate 9, instead of h y , a new state 
£ y , defined as 

& = + 5o >l~h y + 4wX (25) 

is regulated. Thereby, the uncontrollable mode at s = ±ju) 0 and s = ±2 ju 0 are 
embedded in £ y . In a similar way, for regulating h y , a new state defined as 

f* = 9 (V,) + 5 + 4u 4 0 9 (26) 
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is regulated. The result is that h y or 9 become harmonic function with angular 
rates and 2u> 0 in steady state. In this paper, only the design for regulating 0 is 
considered, since the developement of the control law for regulating h y is similar. 
From (22) £ y satisfies 


lyn u y 


(27) 


By defining a state vector x y as 

o o 0( //; ) ^( /v ) ew) 0( v n £ y ^ y 

(21) and (27) can be represented as 


Xy = 


Xy = (Aoy + AA y )x y + (B 0y + A5y)u y (2S) 

where 



T 


Mi 

II 

0 0 

0 

0 0 

-1 0 

lyn 

A B y = 

0 0 

0 

0 0 

< 

i 

0 0 

kyl 15&y n 

4, 

II 

CM 

3&yn — 

5. 
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Note that the pitch angle, 9 is not included in the state vector x y . The variations, 
AA y and A B y , can be decomposed as 


Ay4 y = DyL a y{n&ky)Ey 

A By = FyL b y(Afy)Gy 


where 


Dy = 


Ey = 


0 0 0 0 0 1 0 0 


12 ^ 0 15 s * 0 3 ^ 0 0 0 

Tl 71 Tl 


Fy= - Dy , Gy = [1] 

L a y = [flAty], L\yy = [A/y]. 


The parameter n in above equation denotes the weighting between A k y and A/ y . 

The control law can be obtained by identifying (nAAr y , A/ y ] r as e in ( 2 ), dropping 
the subscript 'y\ and using (9) and ( 11 ) with appropriate choices of p , 7 , n, C and 
C\. Then, the control law u y is represented in the form of 

Uy = KyXy (29) 

where K y is control gain matrix. From the definitions of u y and (29) becomes 

+ 5u try 11 ' + 4u > 4 0 Ty) = K\0 + K\ 9 (30) 

*yn 

+ K% 9 (III) + K » 9 {IV) + Ki 9 {V) + I<1 6 (VI) 

+ + 3uH, + *u%) 

+ 

where Kf denotes the i-th element of the gain matrix K y . The above form is not 
realizable since it needs derivatives of 9 and h v . Fig. 1 describes equation (30). Define 
a new variable x y as 

Xy -j -T y - Ki 9 — K\ 9 - iq f hydt- Ki h y . (31) 

*yn J 
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Then, (30) becomes 


v m + 5 . 2 V (M) , 4 , ,4 * 

= (A7 - 5^A' 6 v )0< n/ ) + (A'S - 

+ (A'jf - 4^A|)<? + - 4w<A'|)0. (32) 

The above equation can be implemented by using the canonical realizations such as 
the controller canonical realization, the observer canonical realization, and the parallel 
canonical realization 20 . In this paper, the parallel canonical realization is adopted. 
Introduce variables £ y and rj y such that 

Cy + Cy = 0 (33a) 

Tj y +4u 2 0 r, v = d. (33b) 

Then, Xy can be represented in terms of £ and tj as 

Xy = Ay(y + ByCy+C y T) v + V y ri v (34) 

where 

B, = 

C, = - J (1 eulKS - 4 Kl + 

v, = - i (vjJ'Kl - 4 k; + a ,;*KZ) . 

Combining (31) and (34) yields an implementable form for the control torque, T y as 

t, = it; e + [<i e + KiJ h„ n + k\ h, ( 35 ) 

+ R\ f, + K\ 4 + Kl n , + Kl i, 

where 

K\ = I yn I<l, I < y 2 = I yn I<l I<l = I y J< y 7 , Kl = I yn I<l 
I<i = I yn A y . I< y e=I yn B y , K y 7 =I yn C y , f<l = I yn V y . 

The control torque discribed by (35) forms a dynamical feedback control law which 
has the same form as in Ref. 2. Fig. 2 describes the control law (35). The integral 
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feedback in (35) is expected to reject the constant input disturbance as in the classical 
control theory 21 . Equation (33) represents the cyclic disturbance rejection filter for 
attitude hold in pitch-axis. The initial states of the integrator in (35) and the cyclic 
disturbance rejection filter are the designer’s choice. 

For the controller design, p , 7, C\, and C are chosen as 

p = 0.81, 7 = 0.2, n = 5, C x = 0 

C = diag ( 3.9a;® 3.9 u 5 0 3.9 u* 3.9u 3.9u; 2 3.9u; 0 ^ ) 

and the minimal nonnegative definite solution to the corresponding ARE is taken. 
Table 2 shows the controller gain matrix K y and the closed loop eigenvalues. A 
stable region for the system parameters of the game theoretic design shown in Fig. 3 
is obtained by applying the Routh-Hurwitz criterion to the closed loop system for the 
given control law. MATHEMATICAL* software is used to check the Routh-Hurwitz 
criterion. Stability margins in some specific direction are listed in Table 3. The stable 
region of the game theoretic design and the LQR. design in Ref. 2 are compared in 
Fig. 3. The bound k y = 0 comes from the open loop characteristic. Fig. 3 shows 
that the game theoretic design improves the stability robustness with respect to the 
parameter variations. A simulation is performed with parameter set considered in 
Ref. 2 

A d ly = 2 ft-lb, A\ y = 0.5 ft-lb, B d = 4 ft-lb, 

<^i y = 0 deg, <^(0) = 1 deg, <^(0) = 0.001 deg/sec, 
other initial conditions = 0. 

Fig. 4 shows the time responses for the nominal system and a perturbed system with 
5 = 60% in Ax variation denoted by solid line and dotted line, respectively. As 
expected, the attitude approachs TEA attitude, —7.6 degree for nominal system, and 
while the CMG momentum oscillates with zero mean value at steady state. 

C. Roll- Yaw Control 

The controller for the roll-yaw axes can be developed in similar way to the pitch 
axis. 
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Define 


e r = <t> - <p c , e 2 = t!> - rp c 

where <f> c and rp c are the command roll and yaw attitude, respectively, and are assumed 
constant. Representing (14a) and (14c) in terms of e x and e., and differentiating the 
resulting equations along with (16a) and (16c) yield 


where f x 


4'" 1 = - (5 + ) - 4(1 + 5/t, W4'"> (36a) 

- 16fc,uife r + (1 - k r )u 0 

■ («<’"> + ^y; v) + 4u*.) - / A 

e^"! - (i, - 5VM 1 ') + (5/t, - 4)^e('"! (36b) 

+ 4^u;^e r — (1 + k z )u 0 

W' = -S^.h^-^h. (36c) 

+ + 5w'/i'"'! + 4^4.) + /„ Ut 

4' 1 "' = - 5 vlhW - 4 u,fh. (36d) 

- + 5 «;»<'"> + 4u,&) + /„u. 


/* — ^jff , and u r and u* are new control variables defined 
^ = ±(^+50,^"'!+^), i = x,*. 

Jin ' 


as 


The system (36) contains uncontrollable modes at 5 = O(double pole), s = 
±7^o(double pole), and s = ±2ju; 0 (double pole) which arise from the external distur- 
bance torque. This means that the external constant disturbance torque and cyclic 
disturbance torque can be rejected in only two of the four states e xy e zy h xy and h z . 
In a similar way to the pitch control, these uncontrollable modes can be removed by 
changing the regulated variables. Tables 4 and 5 show the combination of two states 
in which the constant disturbance and the cyclic disturbance are rejected, respec- 
tively. As shown in Tables 4 and 5, an uncontrollable mode still exists in some of 
the outputs e r . e-. h xy h z . Note that it is always e x that does not reject the cyclic 
disturbance. However, contrary to the pitch channel where bias in pitch angle can 
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not be regulated, the bias in both yaw and roll angle can be rejected, leaving only 
an oscillation in the roll angle. In this paper, only case I for constant disturbance 
rejection and case 6 for cyclic disturbance rejection are considered. 

In order to reject the constant disturbance torque in the attitude channels, the 
uncontrollable double poles at s = 0 are embedded in the CMG-momentum channels 
by regulating h x and h z instead of h x and h z . Similarly, the uncontrollable double 
poles at s = ±ju> 0 and s = ±2ju 0 are embedded in the roll-attitude and yaw CMG- 
momentum channels to reject the cyclic disturbance torque in yaw-attitude and roll 
CMG-momentum channels. By defining ^ and as 


(36) becomes 


& = 4 /V) + ou>le x + 4u*e x 

(, = hW+5ulhW+4u}k, 


$ tT) = -4*^& + (l -k x )u a 

•(e ™ + ^leW+^%)-f xUx 
4 V//) = (k, - 5 Vjei v ) + (5 k z - 4K<ei //7) 

+ - (1 -(- k z )uj 0 ^ - f z u z 

hi Vn = ~ 5o > 2 0 hW _ 4 u>% + u 0 £ z + I xn u x 

6 = -u 0 {h ( P +5bjlhlJ ir >+4ulh x ) + I zn u z . 

By defining a state vector x as 

u i . «. t. i. 4" n 

4 V) h, h, /.<'"> h<-r f. r. 

(37) can be rewritten as a state-space representation of the form 


(37a) 

(37b) 

(37c) 

(37d) 


x = Ax -f B 


Ux 

U : 


(38) 
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where A = 


A a Oiqx6 

^6x 10 A m 


01 00000000 
00 10000000 

0 k tl 0 0 0 k l2 0 k s 3 0 k lA 

00001 00000 

00000 1 0000 
000000 1 000 
00000001 00 

00000000 1 0 

00 000 000 0 1 

o o k zl 0 k z2 0 k z3 0 k z4 0 


A-m — 


0 10 0 

0 0 10 

0 0 0 1 

0 0 0 0 

0 — 0 — 5u >1 

-4 ul 0 -5u >1 0 


0 0 

0 0 

0 0 

1 0 

0 u> 0 

—oj 0 0 


B = 

02x2 

-/« 

02x6 

0 

02X4 

0 

4 n 


. 

0 




hn 

0 


kx 1 - - 4 k x u 2 oy k x2 = 4(1 - k x )u;l k xZ = 5(1 - k x )uj] 
^•r4 = (l k z \ = — (1 + k z )uj 0j k z2 ~ 4 k z LJ^ 

k z 3 — {5k z — 1 )<*;*, k Z 4 = (k z — 


In above system four system parameters, k xy k zy f xy f z are included. However, one of 
them are represented in terms of others. f z are represented in terms of k Xy k z and f x 
as 

/, n 1 + k z 

* z 7 1 _ lc * z ' 

1 rn 1 A -x 
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For the smal] variations of k x , k y , and / z , the variation of /*, A/., is approximated 


as 


AA ~ AA: X 4* &2^k z + ^ 3 A/ x 

where 

_ [*n (Id" ^n) __ An 1 An 1 + i':n 

1 Mi-W’’ ' Cj- 7»T=iT 

Then, the variation of the system and input matrices, AA and AB, can be de- 
composed as 


= D L a (niAk x , n 2 Ak x ) E , 


AB = F L b {n x Ak x ,n 2 Ak z ,Af x )G 


where L a = diagfo Ak x , n 2 Ak z ), 


Lb = diag(n! AA: Z , n 2 Ak z , Af x ), 


D = 


J 2x2 


= (fe) 


F = 


•»4x2 


1 

0 

0 

0 

1 

0 

0 

0 


G = 


-»2x6 


0 

1 


nj 
0 - 


o 4 


X6 


nj 

0 

1 

1 

1 


-i T 


10 0 0 

0 51 * 

n i "2 J 


02x6 


0 0 0 


n 1 


0 4 s 


"2 


nj 


'4x6 


"1 

0 


n 2 


0 


-*2x6 


The parameters nj and n 2 represent the weightings among the three system parame- 
ter. By the choice of n x and n 2 with the ratio of reciprocal of directional derivatives 
of k x , k j and f x with respect to a particular inertia variation direction, the inertia 
variations listed in Table 1 can be assumed in the design process. The directions 
that are preferable for design are the inertia variation that can be made large before 
reaching a physical constraint. 
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The control law can be obtained by using (9) and (11) with appropriate choices 
of p, 7 , n, . n 2 , C and C\. Then, the control law u y is represented in the form of 




= Kx 


(39) 


where K is a control gain matrix. From the definitions of u x , u y , and (39) 
becomes 

T“ ( TP + 5w*T/ //7) + 4u*Ti) = I<\ U + I<\ U (40) 

* t n 

+ I<' 3 U + I<\ e z + I<\ k z + K\ l z + K\ e[ UI ^ 

+ I<i e { / V) + I<' 9 e + K\ 0 ei VI > + K' n h x + K' n h x 

i = x, z where Kj and Kj denote the j-th element of first row and second row of the 
gain matrix K, respectively. The above form can be changed into the realizable form 
in a similar way in the pitch control. Define a new variable x% as 

Xi = j—Ti - I<\ f e z dt - K\t x - K' 3 e x (41) 

An J 

- Ki e z - K[ 0 e z -K[h x - I<[ 6 h z , 
i = x,z. Then, from the definition of ^ and £ z (40) becomes 

x! V) + 5ufo!'" ) + 4u><*, = K\e, + (K\ - 4o>jKj)«, 

+ ( Kj - 4w;K|o)e, + (K‘ - SujKj) 4'"> (42) 

+ ( K‘, - 5 ulKiJeW + (Kj, - 4 u;k;,)A. 

+ K‘ n h .. + (Kj, - SvlK-M '" > + KUlftl' 1 '', 


i = x,z. Introduce variables ( x , ( z , rj x , and r) z such that 

Cx "b Cx = 

7x + 4 T7 r = 

(« + Cr = 0 - 0c 

+ 4 u;* 77* = rp - 0 C . 


(43a) 

(43b) 

(43c) 

(43d) 
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Then, y, can be represented in terms of £ and 77 as 


Xi = A (-* + Bi Cz + c, TJ : + V, tj 2 (44 

+ Si ( r + -Tj Cl + TJz + H, 7 X . * = X , 5 

where 

A= |(«X- re, + «„-’*;) 

= 2 (^cAlO — 2 ^6 ~ w 0 4 ^) 

ft = - 5 (lfejAj - 4 *r; + .,;%•) 

A = - 72 ( Mu2 ° K io - ‘6^1 + 4 w; 2 Ki-^lC t ) 

ft = jG’A-'s-A-ja + ^ir;,) 

*= - j (At — C’*!,) 

ft = - J (lSwjAj, - 4^3 + 

H, = i (iK\, -^K\,) . 

Combining (41) and (44) yields an implementable form for the control torque, T, as 

T, = /q J ( 4 > - <f> c ) dt + K^(<f> - 4 > c ) + K' 3 j> ( 45 ) 

+ K\ J {'P - 'Pc) dt + K' 5 (xl> - 0 C ) + 

+ /?• h x + K' a h z + /?’ Cx + Ki 0 G + A'li 17, 


+ K \1 '/x + ^13 Cx + K\ A Cx + ^[5 T)x + K\ 6 fjx, 


i = x,z where 


Ki = I,J<U Ki = I,„Ki, = 

K\ = Ii«Kl Ki = I in K[ 0 , K ' 7 = I in K[ 5 , Ki=I in K\ 9 , 


A 9 An* 4 i) K x q Ii n Biy /fj] — ^ ln O , , 


— An Ci, -^12 — An T^i , 


A 13 AnG , K \ 4 — Iin^Fii K\$ — linGii Afjg— An?A- 

The control torque described by (45) forms a dynamical feedback control law which 
consists of a conventional PID control and cyclic disturbance rejection filters. The 
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integral feedbacks in (45) are expected to reject the constant input disturbance in 
attitude. Equation (43) represents the cyclic disturbance rejection filter for the yaw 
attitude and the roll CMG momentum. The initial states of the integrators in (45) 
and the cyclic disturbance rejection filter are the designer’s choice. 

For the roll-yaw channel controller design, p , 7, n u n 2 , C,, and C are chosen as 


p = 0.095, 7 = 0.172, n, = n 2 = 5 

“1 0iqx6 


c = 


, Cl = 0 


2x2 


06x10 ^2 

fti = 4.6 • diag( 2.4^ o.lu/* 0.1u; o 2.4u;J u 6 a 

) 

n ! = 1.5-diag( 7“^ £*£ 

£<■* £•* ft**) 


and the minimal nonnegative definite solution to the corresponding ARE are taken. 
The controller gain matrix K and the closed-loop eigenvalues for roll-yaw channel 
are shown in Table 6. The largest closed-loop eigenvalues are seen to remain close 
to the orbital frequency. The stability margins in some specific variations are listed 
in Table 3. For all type of variations listed in Table 3 except Aj, A 3 , and A 6 , the 
designed controller stabilizes the system far beyond the physical limit which means 
that good performance robustness is achieved for these directional variations. For 
A 3 , and Ag, 62% stability margin is achieved. A simulation is performed with 
parameter set considered in Ref. 2 


A\ x = A d u = lft-lb, A\ x = A\ z = 0.5 ft-Ib, 

Bi = B d = 1 ft-lb, v? lr = ip 2x = <p lx = ip 2x = 0, 

0(0) = 0(0) = 1 deg, <£(0) = 0(0) = 0.001 deg/sec, 

0c = 0c = 0 deg, other initial conditions = 0. 

Fig. 5 shows the time responses for the nominal system and a perturbed system with 
6 = 60% in Ai variation denoted by solid line and dotted line, respectively. With 
no noticable performance degradation, the system appears to have good performance 
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robustness. The constant disturbance torques are rejected in roll-yaw attitude chan- 
nels while the cyclic disturbance torques are rejected in roll-CMG and yaw attitude 
channels. The CMG momentum in the roll channel approachs to a constant value 
while the CMG momentum in the yaw channel oscillates around a constant value. 
The biased CMG momentum in steady-state can be changed by changing the com- 
mand attitudes, (f> c and 0 C . The CMG momentum in roll-yaw channel is unbiased 
when the command attitudes are set to the torque equilibrium attitude(TEA). 


IV. Conclusions 

The game theoretic controller is applied to momentum management and atti- 
tude control of the space station in the presence of uncertainty in the moments of 
inertia. The game theoretic controller has been developed for an uncertain linear 
time- invariant system by representing the uncertain dynamic system as an internal 
feedback loop and considering the the input and the fictitious disturbance caused by 
parameter uncertainty as two noncooperative players. It was shown that this con- 
troller stabilizes the system for the prescribed parameter uncertainty bounds. Inclu- 
sion of the external disturbance torque to the design procedure results in a dynamical 
feed back controller which consists of conventional PID control and the cyclic distur- 
bance rejection filter. This shows the state space formulation for design provides a 

v — --- ... 

proper mechanization for handling the external disturbance. It was shown that the 
game theoretic design achieves a stability robustness with respect to inertia varia- 
tions without sacrificing performance robustness, and without increasing the system 
bandwidth. 
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Table 1. Variation type and physical limit of variation due to the triangular 

inequality 

Variation Type 


(A; — [ AI X AI y A I z ]) Physical Bound of 8 


A! = <5[ Ixn Iyn /,„ ] 

-100.0% 

< s < 

oo 

A 2 = <5[-/ xn I yn I zn ] 

-15.9% 

< 8 < 

2.6% 

a 3 = <[/„ /„] 

-81.9% 

<8 < 

13.1% 

A 4 = (S[ 7 rn I yn — I zn ] 

-2.1% 

<6 < 

19.5% 

A s = i[o 

-27.5% 

<8 < 

5.3% 

A 6 = 6[/ rn 0 I zn ) 

-90.0% 

< 8 < 

30.3% 

A 7 = <$[ I xn I yn 0 ] 

-4.1% 

< 8 < 

48.4% 

As = <5[ 0 — Iyn hn } 

-40.0% 

< 8 < 

3.6% 

A 9 = 6[ —I xn 0 I zn ] 

-17.5% 

< 8 < 

2.3% 

Aio = 5[— I xn I yn 0 ] 

-31.3% 

< 8 < 

6.4% 
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Table 2: Controller gains and closed loop eigenvalues for pitch channel 


i Ry 

unit 


1 3.9142E+2 

ft-lb/rad 

m 

2 1.7736E+5 

ft-lb- sec/rad 


3 1.9885E— 6 

ft-lb/ft-lb-sec 2 

m 

4 6.6961E— 3 

ft-lb/ft-lb-sec 


5 2.7357E— 6 

ft-lb- rad/sec 2 

m 

6 4.9017E— 2 

ft-lb-rad/sec 


7 — 2.2665E— 4 

ft-lb-rad/sec 2 

■ 

8 2.0948E— 1 

ft-lb-rad/sec 


Closed loop eigenvalues 

- 

-4.77, -1.52, 

-0.55 ±0.42; 

a 

—0.13 ± 1.00; 

-0.59 ± 1.99; 

* 


Table 3: Stability margin of variation 


Variation 

Lower & upper margin of 6 


Type 

Pitch 

Roll- Yaw 


Ai 

-99% < S < 82% 

-90% <8 < 75% 


a 2 

— 7%° < 8 < 19% 

-90% < <5 < 99% 

m 

A3 

-99% <8 < 99% 

-62% <8 < 56% 


a 4 

-22% <8 < 7%“ 

-73% <8 < 63% 

§ 

A s 

-14%° <8 < 33% 

-99% <8 < 99% 

m 

^6 

-99% < <5 < 99% 

-73% < 8 < 66% 

m 

A 7 

-51% <6 < 16%“ 

-99% <8 < 99% 


a 8 

-14%“ <8 < 43% 

-62% <8 < 66% 


Ag 

— 7%“ <8 < 20% 

-74% < 8 < 89% 

W 

A]n 

—16%“ < <5 < 37% 

-99% <8 < 99% 

it 


“This bound comes from the open-loop characteristic in pitch axis. H ! 
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Table 4: Rejection of the constant disturbance torque for roll-yaw axes 


Case 

States 

Uncontrollable mode 
in resulting system 

1 

^x ? € 2 

none 

2 

h Xy h z 

none 

3 

e X J h x 

none 

4 

e z > h z 

none 

5 

J h z 

5 = 0 

6 


5 = 0 


Table 5: Rejection of the cyclic disturbance torque for roll-yaw axes 


Case 

States 

Uncontrollable mode 
in resulting system 

1 

^z 

5 = ±ju 0 

2 

^x ? ^z 

none 

3 

e X ) ^x 

5 = ±ju; 0 

4 

^Z 5 h z 

5 — 

5 

^X ) ^Z 

5 = 

6 


none 
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Table 6: Controller gains and closed loop eigenvalues in roll-yaw channel 


i 

I<? 

i<! 

unit 

1 

i 

8.0779E— 2 

6.9890E-2 

ft-lb/rad-sec 


2 

9.2311E+2 

7.5497E+2 

ft-lb/rad 

m 

3 

4.1496E+5 

2.1476E+5 

ft-lb-sec/rad 


4 

— 1.6675E— 2 

1.2840E-2 

ft-lb/rad-sec 

m 

5 

— 1.7214E+2 

2.1965E+2 

ft-lb/rad 


6 

1.6537E+5 

3.6144E+5 

ft-lb-sec/rad 

m 

7 

2.795SE— 3 

2.1891E— 3 

ft- lb/ft- lb-sec 


8 

1.250SE— 3 

1.6052E— 3 

ft-lb/ft-lb-sec 

« 

9 

— 1.0291E— 4 

— 1.3627E— 4 

ft-lb-rad/sec 2 


10 

-6.9590E-2 

-7.9080E— 2 

ft-lb-rad/sec 

m 

11 

2.4992E— 4 

— 3.3459E— 4 

ft-lb-rad/sec 2 


12 

5.8883E-2 

— 2.4604E— 2 

ft-lb-rad/sec 


13 

— 1.1076E-10 

-1.6908E-10 

ft-lb/ft-lb-sec 3 

fI 

14 

3.0504E— 7 

1.4664E— 7 

ft-lb/ft-lb-sec 2 

• 

15 

16 

-1.0124E— 9 
— 4.6107E— 7 

— 1.3267E— 9 
— 4.3079E— 7 

ft-lb/ft-lb-sec 3 

ft-lb/ft-lb-sec 2 



Closed loop eigenvalues in roll-yaw channel 



-2.62, -1.75, 

-1.05 ±0.06; 



-0.10 ± 0.98; 

-0.28 ± 1.09; 

m 


-0.21 ± 0.06; 

-0.23 ± 0.90; 

EH 


-0.10 ± 1.97; 

-0.38 ± 2.05; 
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